Solitons in nonlocal nonlinear media: exact results 
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We investigate the propagation of one-dimensional bright and dark spatial solitons in a nonlocal 
Kerr-like media, in which the nonlocality is of general form. We find an exact analytical solution 
to the nonlinear propagation equation in the case of weak nonlocality. We study the properties of 
these solitons and show their stability. 



I. INTRODUCTION 

Let us consider a phenomenological model of nonlocal 
nonlinear Kerr type media, in which the refractive index 
change An induced by a beam with intensity I(x, z) can 
be represented in general form as 

/oo 
R(x' -x)I(x',z)dx', (1) 
-oo 

where the positive (negative) sign corresponds to a focus- 
ing (defocusing) nonlinearity and x and z denote trans- 
verse and propagation coordinates, respectively. The 
real, localized, and symmetric function R(x) is the re- 
sponse function of the nonlocal medium, whose width 
determines the degree of nonlocality. For a singular re- 
sponse, R(x) = S(x), the refractive index change becomes 
a local function of the light intensity, An(I) = ±I(x, z), 
i.e. the refractive index change at a given point is solely 
determined by the light intensity at that very point. 
With increasing width of R(x) the light intensity in the 
vicinity of the point x also contributes to the index 
change at that point. In the limit of a highly nonlo- 
cal response Snyder and Mitchell showed that the beam 
evolution was described by the simple equation for a lin- 
ear harmonic oscillator |l| . The influence of nonlocality 
of the nonlinear response on the dynamics of beams was 
illustrated for the special logarithmic nonlinearity, which 
allows exact analytical treatment [^) . 

While Eq. @) is a phenomenological model, it never- 
theless describes several real physical situations. Possible 
physical mechanisms responsible for this type of nonlin- 
ear response includes various transport effects, such as 
heat conduction in materials with thermal nonlinearity 
|^-^| , diffusion of molecules or atoms accompanying non- 
linear light propagation in atomic vapours Q, and drift 
and/or diffusion of photoexcited charges in photorefrac- 
tive materials |7],|8| . A highly nonlocal nonlinearity of the 
form (|l|) has also been identified in plasmas PHl^] and 



it appears as a result of many body interaction processes 
in the description of Bose-Einstein condensates |l3| . 

Even though it is quite apparent in some physical situ- 
ations that the nonlinear response in general is nonlocal 
(as in the case of thermal lensing), the nonlocal contri- 
bution to the refractive index change was often neglected 
p^ , p^[ . This is justified if the spatial scale of the beam 
is large compared to the characteristic response length of 
the medium (given by the width of the response func- 
tion). However, for very narrow beams the nonlocal- 
ity can be of crucial importance and has to be taken 
into account. For instance, it has been shown theoreti- 
cally that a weak nonlocal contribution arrests collapse 
(catastrophic self-focusing) of high power optical beams 
in a self-focusing medium and leads to the formation of 
stable 2D (diffracting in two transverse dimensions) soli- 
tons |ll|]L^ JH||l7| . On the other hand, a purely nonlocal 
nonlinearity leads to formation of so-called cusp solitons, 
which, however, are unstable 

Some of the consequences of nonlocality in the nonlin- 
ear response have been observed experimentally. Suter 
and Blasberg reported stabilization of 2D solitary beams 
in atomic vapors due to atomic diffusion, which carries 
excitation away from the interaction region @. Also, the 
discrepancy between the theoretical model of dark soli- 
tons and that observed experimentally in a medium with 
thermal nonlinearity has been associated with nonlocal- 
ity of the nonlinearity [ ^4|JT^1 . 

Here, we investigate the propagation of ID beams in 
nonlinear media having a weakly nonlocal response of 
the general form (Q). Our goal is to find exact analyti- 
cal solutions for bright and dark spatial solitons and use 
them to determine soliton properties, such as existence 
and stability. We start with the paraxial wave equation 
describing propagation of a ID beam with envelope func- 
tion ip = ip{x, z) and intensity / = I{x, z) — {^(x, z)\ 2 , 

id z ^> + hefal) + An(I)i> = 0. (2) 
When the nonlocality is weak, i.e. when the response 
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function R(x) is narrow compared to the extent of the 
beam, we can expand I(x', z) around the point x' = x to 
obtain 

An(I)=±(I + 1 d 2 I), (3) 
where the nonlocality parameter 7 > is given by 
1 1"°° 

7 = - / R(x)x 2 dx, (4) 
1 J -00 

and where we have assumed that the response function is 
normalized, f_ R(x)dx — 1. Note that for R(x) = S(x), 
7 = and Eq.°([|) describes the local Kerr nonlinearity. 
For weakly nonlocal meda 7 <C 1 is a small parameter. 
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FIG. 1. Intensity profile I(x') — I(x', z) and response func- 
tion R(x' — x) in the weakly nonlocal limit. 

Substituting An(I), given by Eq. (^), into Eq. (||) gives 
the modified nonlinear Schrodinger equation 

id^ + ±Efcl> ± (|VI 2 + id 2 x \M 2 )i> = 0. (5) 

The weak nonlocality appears thus as a perturbation to 
the local nonlinear refractive index change. For a single 
peak beam in a self-focusing medium this perturbation 
is of negative sign in the central part of the beam, where 
it serves to decrease the refractive index change. Hence, 
even for very narrow and sharp intensity distributions, 
the resulting self-induced waveguide will be wide and a 
smooth function of the transverse coordinates. In some 
sense, this is similar to saturation of the nonlinearity. 
One may therefore expect that certain features of soli- 
tons of Eq. (H) will be reminiscent of those exhibited by 
solitons in saturable nonlinear media, ft transpires, how- 
ever, that nonlocality also results in new effects, espe- 
cially for self-defocusing nonlinearities. In the following 
we consider separately the case of self-focusing and self- 
defocusing nonlinearities. 

II. BRIGHT SOLITONS 

For a self-focusing nonlinearity the sign of the refrac- 
tive index change is positive. We search for a stationary 
bright soliton solution to Eq. (||) of the form 



ip(x, z) = u(x) exp(irz), (6) 

where the profile u = u(x) is real, symmetric, and expo- 
nentially localized and the propagation constant T > is 
positive. For this solution Eq. (||) reduces to 

d 2 x u + 2{u 2 -T)u + 2 1U d 2 x u 2 = 0, (7) 

which can be integrated once to give the equation 

(1 + 4 1U 2 )(d x u) 2 + {u 2 - 2T)u 2 = C, (8) 

where C is an integration constant. For exponentially 
localized solutions C = 0. Using that u(x) has its max- 
imum Mo at the center x — we futher obtain the well- 
known relation between the propagation constant T and 
the amplitude uq 

ul = 2T. (9) 

Equation (^) can therefore be simplified to 

{d x uf = u 2 {u 2 -u 2 )/{l + ± 1U 2 ). (10) 

Interestingly, had the local and nonlocal contributions 
been of opposite signs, i.e. 7 < (as can happen for in- 
tense laser beams in plasmas) , then a solution to Eq. jlfj| ) 
would only exist if the peak intensity po = uj^is less than 
the critical value p CI = (see also |ic||). A final 

integration leads to 

±x = — tanlT 1 f — ) + a/47 tan -1 (V^yo-) ( n ) 

where we have defined the intensity p = u 2 and the nor- 
malized intensity a 2 — (pa — p)/(l + 47/5). This implicit 
relation gives the profile of bright spatial solitons prop- 
agating in weakly nonlocal Kerr-like media. In the local 
limit, 7 = 0, we recover from Eq. (|ll]) the well-known 
profile u(x) — UQsech(uox) of the ID bright soliton in a 
Kerr medium. 
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FIG. 2. Intensity profiles of bright solitons with unit am- 
plitudes po = 1 for different degrees of nonlocality 7. 
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In Fig. |2j we show the intensity profile of the solution 
( |ll| ) for different values of the nonlocality parameter 7. 
Evidently, an increase in the degree of nonlocality re- 
sults in an increase of the soliton width - nonlocality 
smooths out the refractive index profile thereby leading 
to a broadening of the beam. This effect is more clearly 
seen in Fig. 3 where we plot the width of the bright soli- 
tons (FWHM of the intensity profile) versus the degree 
of nonlocality 7 for different peak intensities; the soliton 
width increases monotonically with the degree of nonlo- 
cality. 
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FIG. 3. Soliton width versus degree of nonlocality 7 for 
peak intensities po — l (solid) and po=10 (dashed). 

An important aspect of any family of soliton solutions 
is their stability properties. For single peak solitons sta- 
bility is determined by the dependence of the power 

/oo 
\iP(x,z)\ 2 dx, (12) 
-00 

on the propagation constant T. Solitons are stable if 
dP/dT > and unstable otherwise For the model 

considered here P(T) can be found analytically to 

P = V7 T Q + l±p^ tan-Hv^). ( 13 ) 
V 4 7 

by using Eq. ([To]) . For weakly nonlocal media with 7 <C 1 
the power is approximately given by 

^ = ^o(l + ^7Po~^7 2 Po + ■■■), (14) 

where P = 2^fp^ is the soliton power in the limit of 
a local response with 7=0. The derivative dP/dT can 
easily be found from Eq. ( |l3| ) and it transpires that the 
power is a monotonically increasing function of the prop- 
agation constant (for 7 > 0), indicating that the solitons 
are stable. 



Interestingly again, had 7 been negative, then the 
bright solitons would exist for sufficiently low peak in- 
tensities, po < 1/|47|, but be unstable for po > 0.7/|4'y| , 
for which dP/dT is negative (see again [|10|). 

To demonstrate the stability of the bright solitons for 
7 > we numerically integrated Eq. (||) using the split- 
step fourier method and the exact soliton solution as ini- 
tial condition. In all simulations (for both bright and 
dark solitons) we used a steplength of dz = 0.001 and a 
transverse resolution of da; = 0.05. Results of the numeri- 
cal simulations are shown in Fig. 4 where we demonstrate 
propagation and collision of two bright solitons with unit 
amplitude po = 1. It is evident that the solitons propa- 
gate in a stable manner. The collision, on the other hand, 
is not completely elastic and causes the soliton amplitude 
and width to oscillate slightly. 
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FIG. 4. Collision of unit amplitude bright solitons in a 
weakly nonlocal Kerr-like medium with 7 = 0.1 

The bright solitons of the nonlocal nonlinear 
Schrodinger Eq. (^) were considered by Davydova and 
Fischuk in the context of focusing of upper hybrid waves 
in plasma [Q . In particular, the existence of bright soli- 
ton solutions and their stability was reported by these 
authors. However, the explicit expression (|lTJ) was not 
given and no numerical confirmation was presented. 



III. DARK SOLITONS 

We now consider the impact of weak nonlocality in the 
case of a self-defocusing nonlinearity, which corresponds 
to a negative sign in Eq.(||). We introduce a new spatial 
variable £ = x — Vz, with V being the soliton transverse 
velocity, and look for stationary solutions of the form 
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ip{x, z) = Vp(C) exp(irz + 



(15) 



Integrating Eq. (19) once more we obtain 



where p(C) is the soliton intensity and </>(£) its phase. 
Substitution of this form into Eq. (||) yields 

d c {p(d c <t>-V)} = 0, (16) 
2p(l - 4 1P )d 2 p - {d QP f - V(«9 C 0) 2 + 8Vp 2 d c <t> 
-8p 2 (T + p)=0. (17) 

We are interested in dark solitons, i.e. solutions with an 
intensity dip o n a constant background. Integrating the 
system ( pi|{L7| ) once we obtain that the soliton back- 
ground intensity p = Uq determines the soliton prop- 
agation constant T and that the center intensity p\ = uf 
determines the transverse velocity V, through the rela- 
tions 



v = pi, r = -p 0f 



(18) 



which are exactly the same as those obtained for a purely 
local response. We further find that the soliton intensity 
p(C) is governed by the equation 



(8 c p) 2 = 4(p - Pl )(p Q - p) 2 /(l - 4 7 p). 



(19) 



Obviously, a solution to Eq. ( jr9| ) exists only if the back- 
ground intensity does not exceed a certain critical value, 



Po < Per - 1/|47|, 



(20) 



which coinsides with the critical peak intensity above 
which bright solitons does not exist in a focusing non- 
local medium when 7 < 0, as discussed above. 

It can be shown that the relation ( po| ) represents the 
stability condition for the plane wave solution to equa- 
tion (^|) with a defocusing nonlinearity. For background 
intensities larger than p cr the plane wave solutions be- 
comes modulationally unstable p9|. 
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FIG. 5. Intensity profiles of dark solitons with p Q = 1 and 
p\ = 0.1 for different values of the nonlocality parameter 7. 



±C=^tanh- 1 C£j + v^tan-^v^), (21) 



which is an implicit relation between the soliton inten- 
sity p and the spatial coordinate £, with the normalized 
intensity now given by S 2 (p) = (p — pi)/(l — 47/5) and 
60 = S(po)- For the soliton phase we obtain 



± d> = tan 



— )- ui a/47 tan X (\/^yS), 
u-ij 



(22) 



where we have used the gauge invariance of Eq. (|^) to 
remove a constant phase contribution (phase in the cen- 
ter) . In Fig. H we present examples of the intensity profile 
of the dark solitons for different values of the nonlocality 
parameter 7. These graphs show that increasing nonlo- 
cality decreases the width of the soliton. To better illus- 
trate this effect we plot in Fig. ^| the width of the dark 
solitons (defined as the distance between points where 
the intensity is (pa + pi)/2) versus 7 for several values of 
the soliton contrast a = (p — Pi)/Po- 
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FIG. 6. Width of dark solitons versus degree of nonlocality 
7 for p — 1 and different contrasts a = (p — pi)/p . 

We see that unlike saturation, which tends to increase 
the width of dark solitons, nonlocality has the opposite 
effect. This is due to the fact that nonlocality causes 
the nonlinear index change to advance towards regions 
of lower light intensity. In the case of bright solitons this 
effect leads to soliton expansion. The plots in Fig. || also 
indicate the interesting effect that for strong nonlocal- 
ity, all dark solitons acquire the same width. Actually, 
in this limit, i.e. when 7 approaches the critical value 
7 = l/(4po) the relation (21) can be inverted leading to 
a simple expression for the soliton intensity profile 
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Po[l 
Po, 



acos 2 (woC)]: 



|tioC| < tt/2 
uoCI > t/2- 



(23) 



The solution (£3|) confirms that in this limit all dark soli- 
tons have the same width £ cr = tt/(2uq) independent of 
the contrast. 

With a nontrivial phase structure, these dark soli- 
tons can be represented in the complex plane describ- 
ing the real and imaginary part of the soliton amplitude 
■0(C) = v / /' ex P(*'^)' Any soliton solution is then repre- 
sented by a trajectory in this plane. Figure fj] shows such 
a phase diagram (7 = 0.1) for several values of the soliton 
contrast. The circle represents the background intensity 
po = 1 and nonlocal soliton solutions are plotted using 
solid lines. The dashed line corresponds to the same- 
contrast soliton solution in the purely local case (7 = 0). 
The total phase change across the soliton is given by the 
angle subtended by two lines connecting the origin with 
points of intersection of the soliton curve and the circle. 
One can see directly from Fig. that the total phase 
change across the soliton in a nonlocal medium is always 
less than in the local case. Furthermore, the soliton phase 
profile closely resembles that of a threshold nonlinearity 



1.0 



0.0 



-1.0 



' <s a=l 




\ 


/ 


0.95/ 




/o.84 ^ 














, 1 




1.0 0.0 


1.0 



Re ¥ (Q 

FIG. 7. Real and imaginary part of dark solitons with 
po = 1 in a weakly nonlocal medium with 7 = 0.1. 

The linear stability of dark solitons can be determined 
by considering the conserved quantities of the nonlinear 
Schrodinger Eq. (^) with defocusing nonlinearity, namely 
the rcnormalized Hamiltonian 



H 



[|^| 2 + (|Vf -po) 2 -7(W 2 ) 2 ] dx 



and the renormalized momentum 



Q = 



w )dx - 



(24) 
(25) 



A soliton solution propagating with velocity V corre- 
sponds to an extremum of the Hamiltonian for fixed mo- 
mentum, 8{H — VQ) — 0. It has been shown that the 
stability criterion for dark solitons is determined by the 
dependence of the momentum on the velocity, Q — Q{V) 
|21fl . For dark solitons to be stable the derivative of the 
momentum with respect to the velocity must be positive, 



dQ/dV > 0. 



(26) 



In the case of a nonlocal nonlinearity the expression for 
the Momentum can be evaluated analytically to 



= — 2p a tan 



+ {pa - Pi) 



-[1 + 4 7 (2 Po - p 1 )}J^-tan- 1 (^So) (27) 



The dependence of the momentum on the velocity of dark 
solitons with background intensity po = 1 is shown in 
Fig. |^ for several degrees of nonlocality 7. The mono- 
tonic increase of this function indicates that the dark 
solitons are stable. 
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FIG. 8. Momentum Q of dark solitons with po = 1 versus 
velocity V for different degrees of nonlocality, 7. 

This conclusion is confirmed by direct numerical sim- 
ulations of Eq. (|5|) w ith initial conditions given by the 
exact solution ( |21[]22| ). Fig. ^illustrates propagation and 
collision of two identical dark solitons with po — lOpi = 1 
in a weakly nonlocal medium with 7 = 0.1. While the 
solitons propagate in a stable fashion, their collision is 
clearly inelastic with radiation being emitted from the 
impact area. 
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FIG. 9. Propagation and collision of identical dark solitons 
with po — 10pi = 1 in a weakly nonlocal medium with 7 = 0.1 

In conclusion, we have studied the properties of bright 
and dark one-dimensional spatial solitons in a general 
weakly nonlocal Kerr-like medium, in which the change 
in refractive index due to local and nonlocal contributions 
are of the same sign. We have found an exact analytical 
solution for both bright and dark solitons and used it to 
find the existence and stability regions of both solutions. 

Stable bright soliton solutions were found to exist in 
focusing nonlocal media for all parameter values. Al- 
though they were previously predicted the solution was 
never explicitly written down. The effect of nonlocal- 
ity is in some sense equivalent to that of saturation, - 
to smooth out the index profile and thereby increase the 
soliton width. 

Previously unknown dark soliton solutions were found 
to exist in defocusing nonlocal media for background in- 
tensities below a certain critical value, which correspond 
to the intensity at which the plane wave solutions become 
unstable. We found that nonlocality leads to narrowing 
of the dark solitons, with all dark solitons acquiring the 
same width in the strongly nonlocal limit, independent 
of the soliton contrast. The total phase change across 
the dark soliton was found to be less that in the purely 
local case, with the phase profile resembling that of a 
threshold nonlinearity. 

Both nonlocal bright and dark spatial solitons appear 
to be stable under propagation. Preliminary studies of 
soliton collision revealed their inelastic character in anal- 
ogy to collisions of solitons of nonintegrable models. 

O. Bang acknowledges support from the Danish Tcch- 
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